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ABSTRACT 

We have carried out a comprehensive study of the molecular conditions and star-forming 
activities in dark cloud LI 174 with multi-wavelength data. Mapping observations of LI 174 in 
13 CO 7 = 2-1 and 12 CO J = 3-2 were performed using the KOSMA 3-meter telescope. Six 
molecular cores with masses ranging from 5 to 31 M and sizes ranging from 0.17 to 0.39 pc 
are resolved. Large area ahead of a Herbig Be star, HD 200775, is in expanding and core 1 is 
with collapse signature. Large line widths of 13 CO / = 2 - 1 indicate the ubiquity of turbulent 
motions in this region. Spectra of 12 CO J = 3 - 2 prevalently show conspicuously asymmetric 
double-peaked profiles. In a large area, red-skewed profiles are detected and suggestive of a 
scenario of global expansion. There is a large cavity around the Herbig Be star HD 200775, the 
brightest star in LI 174. The gas around the cavity has been severely compressed by the stellar 
winds from HD 200775. Feedbacks from HD 200775 may have helped form the molecular 
cores around the cavity. Seventeen 2MASS potential young stellar objects were identified 
according to their 2MASS colour indices. The spatial distribution of the these 2MASS sources 
indicates that some of them have a triggered origin. All these suggest that feedbacks from a 
Herbig Ae/Be star may also have the potential to trigger star forming activities. 
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1 INTRODUCTION 

LI 174 is a molecular cloud located in the constellation of 
Cepheus with low- to inte rmediate-mass star forming activities 
l lKun. Kiss. & Balod 12003 ), As a portion of the L1167/L1174 
comple x, LI 174 was reported for the first time as a dark neb- 
ula by iLvndsl dl962b . Investigations of LI 174 have been per- 
formed with multi-wavelengths continuum and emission lines ob- 
servations. The firs t mapping observations with m olecular lines 
were carried out by iMvers. Linke. & Benson! d 19831) in 13 CO J = 
1 - and C 18 J = 1-0. Studies with lines from 12 CO 
and it s isotopic variants have b e en performed by Myers et all 
dl988l).IWu. Zhou. & Evans! Jl992h.lButner. Lada. & Lorenldl995hT 
iBontemps et"al] Jl996h. I Williams etalj dl998b . iBuckle & Fulien 1 
1 20021) . and IWalsh. Myers. & Burtonl d2004l) . Wing emissions of 
12 CO 7=1-0 were detecte d in LI 174 and interpreted as ev- 
idence of outflow activities by IWu. Zhou. & Evani d 1 992b- With 
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observations, e.g. 



1997 1 ) , C 3 H 3 b y iMadden _ . 

19891), HC^N by iFuller&MversId 19931), DCO + a n d H 13 CQ + by 



jTurner. Pirogov. & M inh 
dl989h . CS by IZhou et"al 



observations in NH 3 (J,K) = (1. 1). iGoodman et al.l dl993h de 



tected an averaged velocity gradient of 0.87 ± 0.32 km s 1 
pc~' in LI 174. This region has been studied with lines from 
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iButner. Lada. & Lorenl dl995l) and IWilliams et al.l d 19981) , N 2 H + 
and C? H> by iBenson. Caselli. & Mverj d 19981) . and CH 3 OH, c- 
C 3 H 2 bv lBuckle & Fulled d2002h ." Efforts to search for water masers 
in LI 174 obtained negative results dPersi. Palagi. & Fellil [l994l 
iFuruva et alj2003tlSunada et alj|2007l) . 

Located in LI 174, a well known reflection nebula NGC 
7023 has been extensively studied. NGC 7023 is illuminate d 
by a Herbig Be star HD 200775 dKun, Kiss. & Baloj \200j ). 
There are 3 photodissociation regions (PDRs) in NGC 7023 
located to the east, south and northwest of the exciting star. 
During the past decade, chemical and physical conditions of 
these PDRs have been investigated with observations at wave- 
lengths from mid-infrared to millimeter ( An & Sellgrem 120031: 
Sellgren. Uchida. & Wernerl2007l:lBerne et a l. 2008; Se llgren et all 



Sellgren, Ucnida, & Werner zUU / ; Heme et al. zUUo ; Sellgren et al. 
20iq|Boersma et alfeOloHAbergel _ et al.ll2010l ; Ijoblin et alfcoiol 



Habart et alj|201 1| ; iRosenberg et al.ll20 1 lh . Many of these studies 



were dedicated to spectroscopic observations of PAH emission fea- 
tures in this region. The first confirmedly detected fullerene C60 in 



2 J. H. Yuan et al. 



interstellar space was discovered in NGC 7023 dWerner et al .120041 ; 
ISellgren. Uchida. & Werner! 120071; ISellgren et alj|2010b. With "ob- 
servat ions using SPIRE and PACS onboard Hersc/ie/. lAbergel et al.l 
for the first time detected far-inffared/submm filaments in 
thePDRsofNGC 7023. 

HD 200775 is a massive Herbig Be star wit h a spectral type 
of B3 (±l)e and a total luminosity of 15,000 L Q (H ernandez et al 
2004D . The binarity of this so urce was confirmed bv lMonnier et al 



intensities were corrected to main beam temperature scale using the 



d2006l) . lAlecian et al.l (2008) estimated the masses of the primary 
and the secondary stars to be 10.7 and 9.3 M Q respectively, and the 
age of the primary to be about 10 5 yr. At the near vicinity, a bicon- 
ical cavity has been very likely excava ted by an energetic bipolar 
outflow in an earlier evolutionary stage dFuente et"al .1 19981) . In a re- 
cent work. lokamoto et aL I J2009b presented the detection of a flared 
disk around HD 200775. 

The distance to LI 174 from the Sun has not been 
firmly determined. In the litera ture, several value s have 
been quoted, e.g., 440 pc dWu Zhou. & Evansl 1 19921 ; 
iGoodman etai] 1 1 9931; iBontemps et alj Il996h 429!^ 9 pc 
( Ivan den Ancker. de Winter, & Tiin A Duel Il998l) . and 288 



pc dStraizvs et al.l ll992F In this paper, we will follow 
IWu. Zhou. & Evansl dl992h to adopt the distance of 440 pc. 

In this paper, we will present a mapping study of LI 174 with 
13 CO 7 = 2 - 1 and 12 CO J = 3 - 2. The physical conditions 
and kinematics are profoundly discussed. With combination of mil- 
limeter and infrared data, star forming activities in LI 174 are well 
investigated. This paper is arranged as follows: we present a de- 
scription of the observations and archival data in Section [2] while 
the results and some preliminary analysis are presented in Section 
[3j in Section [4] we try to comprehensively discuss the data; the 
findings of this work are summarized in Section|5] 
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2.2 Archival data 

The Spitzer Gould Belt Legacy Survey 2 (pr ogram ID: 30574; 
lAllen & Spitzer Gould Belt Legacy Teamll2007l) carried out obser- 
vations toward LI 174 and its v icinity with the InfraRed Array 
Camera (IRAC; lFazio et alj2004l) and the Mu ltiband Imaging Pho- 
tometer for Spitz er (MIPS;lRieke et alj|2004h onboard the Spitzer 
Space Telescope dWerner et al.ll2004l) in 2006 November and 2007 
February. These data are publicly available on IRSA 3 , where we re- 
trieved the post-BCD (post Basic Calibrated Data) mosaics. They 
are directly used in the qualitative analysis of the distribution of 
dust and its spatial relation with molecular gas. 

Archived data from the 2MASS Point Source Catalog (PSC) 
were also used in our work. To guarantee the r eliability of the data , 
we adopted the source selection criteria from iLi & Smith! d200 5l) . 
These include: a) each source extracted from the 2MASS PSC must 
have a certain detection in all J,H& Ks wavebands; b) only sources 
with a A"i-band signal-to-noise ratio above 15 are selected. Young 
stellar object candidates were identified according to their 2MASS 
colour indices (e.g., [7 - H], and [H - Ks]. More details, see Sect. 
13.4b . The distribution of these potential YSOs and their relationship 
with other sources were investigated. 

Photometric data of 6 IRAS point sources were retrieved from 
version 2.1 of IRAS Point Source Catalog. And their color indices 
and infrared luminosity are quantitatively investigated. Their rela- 
tions with and effects to the molecular cloud are discussed in fol- 
lowing analysis. 



2 OBSERVATIONS AND DATA ACQUISITION 
2.1 KOSMA observation 

A region of 23' x 15' around LI 174 was mapped in 13 CO 7 = 2-1 
(220.399GHz) and 12 CO 7 = 3 - 2 (345.796GHz) using the 
Kolner Observatorium fur SubMillimeter Astronomie (KOSMA) 
3-m telescope 1 on Gomergrat near Zermatt in Switzerland. The ob- 
servations were carried out with a 1' x 1' grid using the on-the- 
fly (OTF) mode on April 16th 2003. The reference point was at 
(a = 2l''00" ! 22U3,<5 = 68°12'52".9, J2000). 

The 13 CO 7 = 2 - 1 and 12 CO 7 = 3-2 lines 
were simultaneousl y observed with the dual-channel 230/345 
GHz SIS receiver dGrafetal.1 fi998l) . whose noise temperature 
was about 120 K. The on board ac ousto optical spectrometer 
dSchieder. Tolls. & Winnewisserlll989l) had 1501 and 1301 chan- 
nels with bandwidths of 248 and 442 MHz at 230 and 345 GHz 
respectively. This resulted in channel widths of 165 and 340 kHz 
corresponding to velocity resolutions of 0.22 and 0.29 km s -1 at 
230 and 345 GHz. The system temperatures, T„, s , were 186 and 
245 K at 230 and 345 GHz. The beam sizes at 230 and 345 GHz 
were determined to be 130" and 80" with continuum cross scans 
on lupiter. The forward efficiency F e ff was 0.93 during our obser- 
vations. And the main beam efficiencies B e ff were 0.68 and 0.72 at 
230 and 345 GHz. The pointing was frequently checked on lupiter 
and better than 10". 

The data were reduced and visualized with CLASS and GREG 
programs of GILDAS software dPety|2005l) . The baselines were fit- 
ted with one order polynomial and removed for each line. The line 



3 RESULTS 

Figure[T]shows the I3 CO 7 = 2-1 and 12 CO 7 = 3-2 mapping grids 
of LI 174 and its vicinity. 13 CO 7 = 2-1 spectra show mono-peak 
resembling gaussian profiles. In contrast, profiles of 12 CO 7 = 3-2 
show complicated structures. This suggests that 12 CO 7 = 3-2 
transition is optically thick in this region, while 13 CO 7 = 2 - 1 is 
relatively optically thin. Thus, the molecular gas traced with 13 CO 
7 = 2-1 line is less affected by optical depth. We will use this line 
to derive the physical properties of the cloud. 



3.1 Molecular cores 

As shown in Figure[T](a), the intensity of 13 CO line varies through- 
out the region. This suggests a complex distribution of molecular 
gas on the projected plane of sky. There exists a cavity at the mid- 
west. And relatively denser gas is surrounding the cavity. In the 
following analysis, six molecular cores with different physical con- 
ditions are resolved. 



1 The KOSMA telescope has been transferred to Yangbajing, Tibet, China. 

2 http://www.cfa.harvard.edu/gouldbelt/ 

3 NASA/ IPAC Infrared Science Archive (IRSA) is part of the Infrared Pro- 
cessing and Analysis Center (IPAC), which is operated by the Jet Propulsion 
Laboratory, California Institute of Technology, under contract with the Na- 
tional Aeronautics and Space Administration, http://irsa.ipac.caltech.edu/ 
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Figure 1. Mapping grids of LI 174 and its vicinity, (a) For the line of 13 C0 7 = 2-1. The velocity range for each spectrum covers -5 to 10 km s~' , while the 
intensity ranges from -1 to 13 K. (b) For the line of 12 CO 7 = 3-2. The velocity range for each spectrum is as the same as that in panel (a), while the intensity 
ranges from - 1 to 6 K. 
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Figure 2. Velocity-integrated intensity map of 13 CO 7 = 2—1 emission. The velocity interval covers -2 km s to 7 km s . The contour levels range from 
4.8 to 13.1 K km s with a step of 1.5<x (l<x = 0.69 K km s ). The peak positions of the six molecular cores identified in this paper are labeled with crosses. 
The asterisks mark the candidate YSOs from 2MASS PSC, while the filled triangles indicate IRAS point sources in the field of view, whose corresponding 
identities are listed in Table|3] The red star represents the Herbig Be star HD 200775. 



3.1.1 Identification 



3.1.2 LTE parameters 



Figure [2] shows the velocity-integrated intensity map of 13 CO 7 = 
2-1 rotational transition. We, by eye, identified six intensity peaks 
as cores, which are labeled as core 1 to core 6 with ascending right 
ascension order. The coordinates of the cores were converted from 
the offsets of intensity peaks to equatorial system (J2000) and are 
given in Table[T] The uncertainties of them are referred to the point- 
ing error, which is better than 10". The radius of each core was 
determined according to the nearest point with 50% of its peak in- 
tensity and is presented in Table|2] We estimate that more than 80% 
of the masses of the cores are involved in spheres with these radii. 

Core 1 to 5 are distributed surrounding a cavity, which has 
been excavated by an energetic and bipolar outflow in an ear- 
lier evolutionary st age of the massive Herbig Be star HD 200775 
dFuenteetal.ll 19981) . And core 6 is situated at the most east, about 
8' (~ 1 pc) away from HD 200775. Core 1 is associated with IRAS 
20597+6800, while core 3 is coincident with a 2MASS YSO can- 
didate (2MASS J21013520+6810086). No IRAS point sources or 
2MASS YSO candidates are detected at the near vicinities of the 
other four molecular cores. Core 3, core 4 and core 5 show notable 
intensity gradients at the directions towards HD 200775. This could 
be attributed to the feedbacks from the central Herbig Be star. 

The spectra of 13 CO J = 2 - 1 and 12 CO J = 3 - 2 at the 
peak of each core are shown in Figure [3] Rotational line of 13 CO 
7 = 2-1 shows similar mono-peak profiles at all six cores, while 
that of I2 CO 7 = 3-2 differs from each other. For core 1, line of 
12 CO 7 = 3-2 indicates a blue-skewd profile and broad line wings. 
This line shows mono-peak with a blue wing at core 2 and 4, while 
weak blue-skewed profile feature at core 3 and 5. Contrast to core 
1, a red-skewed profile is observed at core 6. 



To investigate the nature of the cores, it is firstly necessary to obtain 
estimates of some physical parameters, e.g., temperatures, masses. 
The derivation of these parameters follows a rotation temperature 
- column density analysis under t he assumed conditio ns of local 
thermodynamic equilibrium (LTE) jGarden et al ]|l99lh . 

One component gaussian fittings to molecular lines at peaks 
of the cores were carried out. In TableQ] the antenna temperatures, 
velocity centroid, and full widths at half-maximum (FWHM) of 
lines of 13 CO 7 = 2 - 1 and 12 CO 7 = 3 - 2 are given. 

According to iGarden et al.l dl99ll) . under conditions of LTE, 
the measured main beam temperature T„,t and the column density 
of 13 CO {Na co ) can be expressed as follows, 



hv 



-[- 



TYia 



K L exp(^)-l exp(^)-l 

3k exp[/iS7(7 + 1)/k7„] 
Sn 3 Bp 2 (7+ 1) 

(T ex + IiB/3k) 



x [1 - exp(-r)]/ (1) 



,-)] I 
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[1 - exp(-hv/KT, 

where T ex is the exciting temperature which we estimate to be about 
20 ± 4 K from the averaged dust temperature of the 6 IRAS sources 
(see Table[3}, T bg = 2.75 K is the temperature of the cosmic back- 
ground radiation, / is the beam-filling factor which is 1 in our ob- 
servations, B and /j are the rotational constant and permanent dipole 
moment of 13 CO, 7 is the rotational quantum number of the lower 
state in the observed rotational transition. 

The combination of equation Q] and [2] would lead to a revised 
expression of Nn co as, 



3k 2 exp[hBJ(J + l)/KT ex ] (T„ + hB/3 K ) 

X [1 - exp(^)] 
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Figure 3. Profiles of I3 CO J = 2 - 1 and 12 CO 7 = 3-2 emission lines at 
the peaks of the six cores. The red dashed line shows the zero level intensity, 
while the green dash-dotted line indicates the system velocity V™ s = 2.58 
km s . 



exp(Av/«r eJ )— 1 exp(/jv/KT/, ? )-l 



rdv 



1 - exp(—r) 
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where t is the optical depth of 13 CO 7 = 2-1. Given the opti- 
cally thin feature of 13 CO 7 = 2 - 1 emission, the term 



l-expi-T) 



is approximate 1. The velocity-integrated intensities fT m bdv for 
the cores are resulted from the gaussian fitting and given in Ta- 
ble^ To obtain the column densities of molecular hydrogen N w ,, 
we adopt a canonical [CO]/[H2] ab undance ratio of ^ 10 4 , as 
measured in nearby molecular clou ds ( Frerking. Langer. & Wilsonl 



ll987l ; |Pineda. Caselli. & Goodmanll2008l) . Given the near distance 
of LI 174, the isoto pe ratio [ 12 C1/[ 13 C1 in the local ISM is applied, 
which is about 77 [ (Wilson &Roodll 19941) . With an assumed error 
of 10 percent for the isotope ratio, uncertainties are counted while 
deriving the N w , . The resulted column densities of 13 CO and H 2 are 
given in column 4 and 5 of Table|2] 

The finally derived mass of each core is reached with an as- 
sumption of a homogenous sphere structure. 



Mr 



-nR'n H n g m{H 2 ) 



(4) 



where /j g = 1.36 is the mean atomic weight of gas, m(Hj) is the 
mass of a hydrogen molecule, and n H = N H J2R is the number 



density of molecular hydrogen. The resulted LTE masses of the six 
cores are presented in column 7 of Table [2] The masses vary from 



3.1.3 Virial masses 

To collapse for a core, the mass should be larger than its vi rial 
mass. We follow iMacLaren. Richardson. & Wolfendalel dl988l) to 
estimate the virial masses of the cores with the assumption of con- 
stant density distributions. 
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where Av is the FWHM of the observed molecular line (see Table 
[]} and R the estimate of radius given in Tabled Column 8 of Table 
[2] shows the resulted virial masses ranging from 171.2 to 463 M Q 
which are evidently larger than the LTE masses. 



3.2 1 2 CO 7 = 3 - 2 mapping 

Figure [4] (a) shows the velocity-integrated intensity map of 12 CO 
7 = 3-2 (contours) overlayed on that of 13 CO 7 = 2-1 (gray). 
The distribution of 12 CO is more extended than that of 13 CO. This is 
because that emissions from 12 CO trace more diffuse and external 
part of molecular clouds. The positions of intensity peaks of 12 CO 
are not consistent with that of cores traced by 13 CO. Self-absorbtion 
of 12 C0 7 = 3-2 at densest region could explain this phenomenon. 
The large column density at the immediate vicinity of cores makes 
12 CO 7 = 3-2 rotational line to be optically thick. This would 
cause deviation in the observed distribution of l2 CO. 

Another interesting feature from the observations is the di- 
versity of profiles of 12 CO 7 = 3-2 line spectra (see Figure Q] 
(b)). Conspicuously asymmetric double-peaked spectral lines are 
detected in large area of this region. Such asymmetric double- 
peaked profiles are always regarded as effective tracers of kinemat- 
ics of molecular clouds. Asymmetrically blue profiles, lines with 
the peaks skewed to the blue sides, are frequently interpreted as ev- 
idence of inward motions . Such inw ard motio ns could be collapse 
dZhou et alJI 19931 : lEvans et alj|2005l) or infall l lWu & Evans]|2003l : 
IWu et alJl2061T By contrast, the red-asymmetry is observed when 
a centrally concentrated system (i.e., with excitation temperature 
decreasing outward) is expanding. 

Figure |4](fo)-(d) show the spectral profiles of 12 CO 7 = 3-2 
in three subregions labeled with red boxes in panel (a). Lines in 
subregion /, corresponding to core 1, show dominant blue profiles, 
which could be due to infall motions (see Sect. 14. Lit . In a rela- 
tively lager area (about 1 pc in diameter), red profiles are detected 
in subregion //. In our further analysis, this would be attributed to 
the global expansion rather than outflow (see Sect. l4.2t . Much more 
complex line profiles are observed in subregion ///, in which blue 
and red profiles show up alternately. The backflow of stellar winds 
and feedbacks from two potential 2MASS YSOs could induce com- 
plexity to the line profiles in subregion /// (see Sect. 14. 2\ . 



3.3 Associated IRAS sources 

There are 6 IRAS point sources detected in the mapping region 
of KOSMA observations. They are marked in Figure|2]and Figure 
[4] with filled triangles and labeled as IRAS 1 to 6 with ascending 
right ascension. Their genuine IRAS names and IRAS photometric 
data from the catalog are given in Table [3] As shown in Figure [2] 
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Table 1. Fitting parameters of cores. 











13 CO (J=2-l) 






CO (J=3-2) 




Label 


R.A. 


DEC. 


n 


vlsr 


FWHM 


n 




FWHM 




(J2000) 


(J2000) 


(K) 


(km/s) 


(km/s) 


(K) 


(km/s) 


(km/s) 


core 1 


21:00:22.13 


68:12:25.9 


2.7395 


2.638(0.014) 


2.234(0.036) 


5.2106 


2.401(0.032) 


3.974(0.083) 


core 2 


21:00:45.84 


68:05:30.7 


2.5911 


2.287(0.014) 


2.126(0.033) 


7.9352 


2.135(0.014) 


2.656(0.035) 


core 3 


21:01:31.11 


68:11:46.3 


4.5975 


2.363(0.014) 


2.116(0.034) 


5.3655 


1.724(0.021) 


3.077(0.046) 


core 4 


21:01:31.76 


68:06:39.1 


4.2942 


2.439(0.009) 


1.843(0.023) 


8.4209 


1.911(0.017) 


3.049(0.040) 


core 5 


21:02:17.13 


68:09:52.9 


5.4667 


2.543(0.010) 


2.197(0.024) 


7.9722 


2.500(0.018) 


3.437(0.042) 


core 6 


21:02:51.30 


68:10:54.7 


3.0554 


2.765(0.015) 


2.385(0.038) 


9.3281 


3.115(0.013) 


3.080(0.041) 



Table 2. Physical Parameters of cores. 



Label 


R 


/ T m bdv 


Nl3 C0 


N„ 2 


n («2) 


M LTE 


M VI > 


Pex/K 




(pc) 


(Kkms -1 ) 


(10 15 cm" 2 ) 


(10 21 cm" 2 ) 


(10 3 cm" 3 ) 


(M ) 


(Ms) 


(10 5 KcrrT 3 ) 


core 1 


0.17 


8.9(0.2) 


4.9(0.2) 


3.8(0.5) 


3.6(1.0) 


4.9(2.2) 


174.3(28.1) 


3.2(2.5) 


core 2 


0.18 


8.0(0.1) 


4.4(0.2) 


3.4(0.5) 


3.0(0.8) 


5.0(2.1) 


171.2(25.6) 


2.5(1.8) 


core 3 


0.20 


14.2(0.2) 


7.8(0.3) 


6.0(0.9) 


4.8(1.2) 


10.9(2.1) 


187.6(26.1) 


3.5(2.9) 


core 4 


0.26 


11.5(0.1) 


6.4(0.3) 


4.8(0.7) 


3.0(0.7) 


15.0(1.8) 


183.4(19.8) 


1.5(0.7) 


core 5 


0.24 


17.5(0.2) 


9.7(0.4) 


7.4(1.0) 


5.0(1.1) 


19.3(2.2) 


247.7(26.9) 


3.6(1.8) 


core 6 


0.39 


10.6(0.1) 


5.9(0.3) 


4.5(0.6) 


1.9(0.4) 


31.0(2.4) 


463.0(40.3) 


1.7(0.6) 
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Figure 4. (a) Velocity-integrated intensity map of 12 CO 7 = 3-2 (contours) overlayed on that of 13 CO 7 = 2-1 (gray). The velocity interval covers -2 km 
s to 7 km s" 1 . The contour levels range from 14.6 to 39.3 K km s with a step of 1.5<x (ltr = 1.83 K km s~'). The red dashed boxes locate 3 subregions 
with complex line profiles. Other symbols are the same as the ones in Figure|2] (b) Line profiles of subregion I labeled in panel (a). The offsets are indicated 
at the up-right corner of each sub-panel, (c) The same as (b), but for subregion II. (d) The same as (b), but for subregion III. 
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Table 3. Properties of associated IRAS sources. 





IRAS Name 


Pl2 

(Jy) 


F73 
(Jy) 


^60 

(Jy) 


F 100 

(Jy) 


f a 






(L Q ) 


t£ 

(K) 


2MASS Association 


IRAS 1 


IRAS 20597+6800 


0.250 


0.823 


4.48 


55.13 


1331 


0.52 


1.25 


7.85 


18.96 


Null 


IRAS 2 


IRAS 20599+6755 


0.677 


0.595 


18.98 


1095.00 


1231 


-0.06 


1.45 


122.94 


14.53 


Null 


IRAS 3 


IRAS 21009+6758 


26.660 


76.770 


637.70 


1095.00 


3333 


0.46 


1.38 


364.98 


31.02 


J21013520+6810086 


IRAS 4 


IRAS 21015+6757 


1.524 


4.252 


48.07 


1095.00 


3311 


0.45 


1.50 


133.86 


16.90 


Null 


IRAS 5 


IRAS 21023+6754 


0.266 


0.387 


0.40 


15.46 


3311 


0.16 


0.18 


2.31 


15.46 


J21025943+6806322 


IRAS 6 


IRAS 21025+6801 


0.357 


0.796 


10.66 


25.17 


3333 


0.35 


1.48 


6.46 


28.13 


Null 



a Flux density quality, 1 for an upper limit, 2 for moderate quality, and 3 for high quality. The four numbers for each source refer to the flux density quality of 
IRAS bands at 12, 25, 60 and 100 jum, respectively. 

^ Dust temperatures were derived from color temperatures T^o/ioo) defined bv lHenning, Pfau.& Altenhofj J 1990h : Tj » 7* C (60/100) = ^ t — ■ 

(3+/3)ln( -jj- >-Hf-r-) 



IRAS 20597+6800 (IRAS 1 in this paper) and IRAS 21009+6758 
(IRAS 3 in this paper) are associated with core 1 and HD 200775 
respectively. 

The infrared flux densi ties of the 6 IRAS sources are derived 
with the following equation JCasolietal.lll986h : 



F(10" 13 WnT 2 ) 



F25 F(,o ^100, 

1.75 x ( + — + — + ), (6) 

0.79 2 3.9 9.9 

where F l2 , F 2 s, F 60 , and F 10 o are the flux densities in Jy at 12 /jm, 
25 yum, 60 /im, and 100 /im respectively. With the distance of 440 
pc, we obtain the infrared luminosity, which is given in column 10 
of Table[3] Color indices of log(^j) and log(^j) are calculated and 
shown in Table [3] either. 

Cross identification with 2MASS point source catalogue was 
performed for each IRAS source. The identification process fol- 
lowed the criteria of: a) an associated 2MASS point source must 
locate in the error ellipse of a IRAS source; b) the reddest one was 
selected if two or more 2MASS sources were found in the very near 
reaches of a IRAS source. None 2MASS associations were identi- 
fied for IRAS 1, IRAS 2, IRAS 4 and IRAS 6. For IRAS 3, only 
one 2MASS source fulfills the first criterion and was determinately 
identified as the near-IR association. Two 2MASS point sources re- 
side in the error ellipse of IRAS 5. The closer one possessing larger 
infrared color indices (i.e, [J - H] and [H - Ks]) was undoubtedly 
identified as the association. The results of the cross identification 
are given in column 12 of Table[3] 



3.4 2MASS YSO candidates 

Highly reliable 2MASS point sources were reexamined accord- 
ing to their colour indices (i.e., [J - H], and [H - Ks]] to pick 
out the potential young stellar objects (YSOs). We dotted all the 
sources on the [J - H] vs. [H - K s ] diagram presented as Fig- 
ure [5] The sources with colour indices fulfilling the relations of 
[J - H] < l.S[H - Ks] - 0.103486 and [H - Ks] > 0.5 are re- 
garded as YSOs dLi& Smith! I2OO5T) . This process resulted in 17 
YSO candidates in the mapping region of KOSMA observations. 
All of them are represented with red "+" in Figure|5] We numbered 
the 17 sources according to their distances to the Herbig Be star 
HD 200775, i.e., 2MASS 1 as the nearest one to HD 200775 and 
2MASS 17 as the farthest. In Figure[2]and[4] they are marked with 
asterisks. Among the 17 YSO candidates , 3 one s were identified to 
be pre-main-sequence stars bv lKun et al.ld2009), and 11 ones were 
identified to be YSO candidates bv lKirk et alj <2009l) . Photometri- 
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Figure 5. [/ - H] vs. [H - K s ] diagram. The 17 YSO candidates are 
marked as red "+" (see Sect. l3.4l for details). Black dots represent 2MASS 
field stars which may be reddened main-sequence dwarfs and giants. 
The solid lines are th e loci of the main-sequence dwarfs and giant stars 
teessell & Bretj|l988h. The a rrow shows a reddening vector of Ay = 5 
mag jRieke & Lebofsky 1985). The dot-dashed line indicates the locus of 
dereddened T Tauri stars (Mever, Calvet. & Hillenbrand 1997). The dashed 
lines, which are drawn parallel to the reddening vector, define the reddening 
band for normal field stars and T Tauri stars. Crosses are over-plotted with 
an interval corresponding to Ay = 5 mag. 



cal data from the 2MASS catalogue and colour indices of all YSO 
candidates are given in column 3 to 8 of Table|4] 



4 DISCUSSIONS 
4.1 Molecular cores 

4.1.1 Evolutionary status 

Diverse profiles of 12 CO 7 = 3-2 are detected in distinct cores. 
This, along with the distribution of the infrared sources (i.e., IRAS 
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Table 4. 2MASS YSO Candidates 



— 


2MASS Name 




m H 




J - H 


H - Ks 


J - Ks 


or* 


D 'Angle* 






(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 




(armin) 


2MASS 1* 


J21013691+6809476 


6.111 


5.465 


4.651 


0.65 


0.81 


1.46 


0.47 


0.00 


2MASS 2 #t 


J21012706+6810381 


12.323 


11.150 


10.417 


1.17 


0.73 


1.91 


0.96 


1.24 


2MASS 3 #t 


J21014672+6808453 


11.792 


10.798 


10.159 


0.99 


0.64 


1.63 


0.95 


1.38 


2MASS 4* 


J21015265+6809520 


14.062 


13.004 


12.356 


1.06 


0.65 


1.71 


1.00 


1.47 


2MASS 5 


J21012774+6808114 


13.999 


13.195 


12.663 


0.80 


0.53 


1.34 


0.94 


1.82 


2MASS 6* 


J21012734+6811383 


15.378 


13.154 


11.806 


2.22 


1.35 


3.57 


0.96 


2.05 


2MASS 7 


J21015475+6806590 


14.116 


13.264 


12.664 


0.85 


0.60 


1.45 


0.87 


3.26 


2MASS 8* 


J21010870+6812525 


14.736 


13.627 


12.902 


1.11 


0.73 


1.83 


0.92 


4.05 


2MASS 9* 


J21021404+6807306 


15.768 


14.656 


13.604 


1.11 


1.05 


2.16 


0.62 


4.14 


2MASS 10* 


J21014391+6814033 


15.036 


13.734 


12.402 


1.30 


1.33 


2.63 


0.57 


4.31 


2MASS 1 1 


J21022228+6812121 


16.434 


15.180 


14.339 


1.25 


0.84 


2.10 


0.89 


4.86 


2MASS 12 


J21023492+6812024 


15.337 


14.377 


13.685 


0.96 


0.69 


1.65 


0.84 


5.84 


2MASS 13 


J21002149+6807452 


17.019 


15.261 


13.934 


1.76 


1.33 


3.09 


0.77 


7.30 


2MASS 14 


J21001891+6811062 


15.386 


14.533 


13.994 


0.85 


0.54 


1.39 


0.98 


7.37 


2MASS 15* 


J21025484+6806210 


15.661 


13.712 


12.261 


1.95 


1.45 


3.40 


0.78 


8.02 


2MASS 16* + 


J21025943+6806322 


13.871 


13.083 


12.357 


0.79 


0.73 


1.51 


0.65 


8.34 


2MASS 17* 


J2 1030756+6808339 


12.435 


11.398 


10.732 


1.04 


0.67 


1.70 


0.95 


8.52 



* Identified to be YSO candidate bv lKirk et al] <2009h. 

' Identified to be pre-main-sequence star bv lKun et alj J2009t) . 

* 2MASS index defined as a = - { grgjjnjq; 1035 ■ where [7 - H] and [H - Ks] are 2MASS color indices. While a = 1, 
a 2MASS point source will be located on the right dashed line in Fig ure \5\ which is the locus dividing red dened normal 
field stars and T Tauri stars l lBessell & Brettlll988l:lMever. Calvet. & Hillenbran3ll997l ; lLl& Smithll2005l) . An younger 
YSO would possess a smaller 2MASS index. 

* Angle distance to HD 200775. 



source and 2MASS YSO candidates), suggests deferent evolution- 
ary status of the cores. 

Core 1 is the smallest one in size. The association with IRAS 
20597+6800 indicates ongoing star formation in core 1. The pro- 
file of 12 CO 7 = 3-2 shows blue-skewed asymmetric feature 
suggestive of infall motions. Bas ed on the simple an alytic model 
of radiative transfer proposed by iMvers et all (U996I) . the inward 
speed is estimated to be about 0.43 km s' 1 . Shown in Figure [6] 
(a) is the position-velocity (P-V) diagram with a slice crossing the 
core 1 from nort h to south. This is a typical P-V diagram show- 
ing blue profile dWu et al1l2007l) . which well resembles an infall 
scenario. Meanwhile, this P-V diagram indicates the detection of 
prominent line wings in core 1. Such feature is always inter preted 
as th e result of outflow activities in star forming regions dLadal 
Il985l) . Point observation of LI 174 (correspon ding to core 1 in this 
study) in 12 CO 7=1-0 was performed by IWu. Zhou. & Evansl 
dl992l) . and the prominent wing emission was interpreted as the ev- 
idence of the existence of an outflow. And in the outflow catalog of 
IWu. Huang. & HeUl99fj) , LI 174 was classified as a source driving 
a bipolar outflow. Shown in Figure \6\(b) is the velocity-integrated 
intensity map of the wing emissions. The red lobe of the outflow 
is unambiguously revealed, while the blue lobe is not clear. Since 
core 1 is located in the large molecular envelope of HD 200775 
whose feedback severely impacts the ambient material, it is diffi- 
cult to separate the blue lobe of the outflow and the surrounding 
gas. 

Core 2 is the lightest one in mass. The 12 CO 7 = 3-2 line of 
core 2 shows mono-peak feature. This core is relatively diffuse. No 
IRAS source or 2MASS YSO candidate is detected in core 2. All 
these are suggestive of an evolutionary stage prior to the beginning 
of star formation. 



Core 3 is associated with a potential YSO 2MASS 
J21012734+6811383 (i.e., 2MASS 6 in this paper). This 2MASS 
source possesses large infrared excess ([7 - H] = 2.22, [H - Ks] = 
1.35) indicative of deeply embedded feature and an early evolution- 
ary stage. The line of I2 CO 7 = 3-2 shows non-gaussian profile 
with the peak skewed to the blue end. This may be attributed to po- 
tential infall motions. As shown in Figure [2] the southern portion 
of this core is compressed by the feedback of HD 200775. 

Core 4 shows a mono-peak profile in 12 CO 7 = 3-2. The 
density is as low as that of core 2. The lack of infrared source in 
this core suggests a similar evolutionary stage as core 2. 

Core 5 is the densest one. The large gradient in Figure ^in- 
dicates that core 5 is heavily compressed by the winds from HD 
200775. In the observations with I2 CO 7 = 3-2, this core is re- 
solved into two cores (see Figure [4). An IRAS source (i.e., IRAS 
21015+6757) is located between core 5a and 5b. Core 5b is associ- 
ated with 2MASS J21021404+6807306 (i.e., 2MASS 9 in this pa- 
per). Compared with 2MASS 6, 2MASS 9 shows smaller infrared 
excess ([7 - H] = 1.1 1, [H - Ks] = 1.05). This indicate that core 
5b could be more evolved than core 3. Shown in Figure [7] (a) is a 
position- velocity diagram along centers of th e core 5a and 5 b. This 
is a typical P-V diagram showing red profile dWu et all2 007). Both 
of the two cores show notable red-skewed asymmetric profiles (see 
Figure[7](Z>)). These features could be ascribed to the global expan- 
sion originated from the feedback of HD 200775 (see Sect. \4.2\ . 

Core 6 is the largest, most massive and diffuse one. The line of 
12 CO 7 = 3 - 2 at the peak of this core shows red-skewed asymmet- 
ric profile. No infrared source is detected in this core. This, com- 
bined with the diffuse feature, suggests that this core is in a very 
early stage. As shown in Figure |4] core 6 resides in the subregion 
// where red-skewed lines dominate. The origin of red-skew profile 
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offset R.A. (arcmin) 

Figure 6. (a) 12 CO J = 3 - 2 position-velocity diagram constructed along 
the vertical dash-dotted line (cross the core 1 from north to south) in Figure 
[4] (a)- The contour levels start from the 20% of the peak intensity with a 
step of 8%. And the peak intensity is 6.18 K. The vertical dash-dotted line 
indicates the system velocity Vjyj = 2.58 km s _I . The two dashed lines 
indicate the beginning of the blue and red wings, (b) Velocity-integrated 
intensity map of wing emissions of 12 CO 7 = 3-2. The velocity interval 
covers -1 to 0.5 km s _I for the blue wing (dashed lines), while 4 to 6 km s -1 
for the red wing (solid lines). The contour levels start from 1.95 K km s -1 
with an increasing step of 0.37 K km s _1 for the blue wing and start from 
2.11 Kkms~' with an increasing step of 0.33 K km s _1 . The cross presents 
the center of core 1, and the filled triangle indicates the location of IRAS 
20597+6800 (IRAS 1 in this paper). 



of 12 CO J = 3 - 2 in core 6 could be due to the global expansion 
(See Sect.l42l. 

4.1.2 Potential collapse 

Compared wit h that in other low- mass star forming regions (e.g., 
case in Taurus, lOnishi et al.ll2002h . the line widths of cores in our 
observations are relatively large, ranging from 1.8 km s" 1 in core 4 
to 2.4 km s _1 in core 6 (FWHM of 13 CO 7 = 2-1 line). In regions 
with temperatures of about 20 K (see Tabled, the isotherm al sound 
speeds would be about 0.27 km s 1 dMcKee & Tar]|2003l) . This is 
much smaller than the observed line widths. Therefore, some other 
non-thermal motions play dominant roles in the broadening of lines 
in this region. Several sources could be responsible for these non- 
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Figure 7. (a) 12 CO 7 = 3-2 position-velocity diagram constructed along 
the vertical dashed line in Figure [4] (a). The contour levels start from the 
15% of the peak intensity with a step of 7.5%. And the peak intensity is 
10.90 K. The vertical dashed line indicates the system velocity V„ s = 2.58 
km s _I . (b) Profiles of 13 CO J = 2 - 1 and 12 CO J = 3 - 2 at the peaks of 
core 5a and 5b. 



thermal line widths. One of them is that the cor es formed in a tur - 
bulent environment and they are still turbulent dSaito et alj|200ot) . 
Meanwhile, some motions after the formation of a central protostar 
would lead to large line widt hs. Such motions could be infall, out- 
flows , winds, and rotations dPavlyuchenkov et al .1120081 : IWu et al.l 
l2005h . Bulk motions due to central protostar may play key roles in 
the case of core 1 and core 3, especially in core 1, where infall of 
material is detected (see Sect. 14.1.11 . Without associated infrared 
sources, the other four cores should be in turbulent conditions. 

From the calculation performed in Sect. 13.11 one can see that 
the derived virial masses largely exceed the LTE masses for all six 
cores. This suggests that there i sn't su ffici ent gravitational en ergy 
to bind the systems. iLee et al.l d201ll) and lWang etail d2008h ob- 
tained similar LTE and virial m asses in IRAS 05 345+31571 and 
MWC 1080 star forming regions. ISaitoetal.ld2006l) also concluded 
that non-turbulent cores have similar virial masses to LTE masses 
but the virial masses are usually larger than the LTE masses for 
turbulent cores. External pressure is needed for a turbulent core to 
maintain a boun d system. 

We follow ISaito et al.l d2006l) to estimate the least external 
pressures for the 13 CO cores to be in bound systems. The virial 
equation will be amended to be as following with external pres- 
sure: 

F = 2U + Q. - 4nR 3 P ex (7) 

where U = jMcr 2 (cr-velocity dispersion) is the kinetic energy, 
Q = - iG 5 "~ is the gravitational energy, and R is the radius of the 
core. The velocity dispersion is deduced from the line width of 
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13 CO using the equation of tr = Av/2 V2/n(2). While F = 0, we 
get the lower limits of the required external pressure. The results 
are then divided by Bolzmann constant k and tabulated in Table[2] 
Notice that all the six cores need external pressures of P ex /k ~ 10 5 
K ci rr 3 to maintain bound systems. 

iMcKee & Tar] 120031) suggested that the pressure at the sur- 
face of a core is related to its surface density: 

P s ,co,e = 0.85 x 10 9 E 2 K cm" 3 (8) 

where S = l.36N H m(H 2 ) is the surface density in g cirr 2 . In our 
case, the resulted surface pressure for all the six cores is about 10 5 
K cnr 3 , which is commensurate with the required external pres- 
sure for maintaining bound systems. In addition, winds from the 
central bright Herbig Be star HD 200775 would provide supple- 
mental pressure to help binding the neighbour cores, i.e., core 3, 
core 4, and core 5. Therefore, we suggest that these turbulent cores 
can be bound with external pressure. 

As discussed above, turbulence plays great role in balanc- 
ing the gravitational force in the cores. However, previous studies 
showed that this support would not exist long enough to prevent 
them collapsing. This could be attributed to the relatively short dis- 
sipation timescale for turbulent motions, which is about ~ 10 5 yr, 
much shorter than the free-fall timescale of ~ 10 6 yr for a core 
with n m ~ 10 3 cnr 3 iMcKee & Ostrikej|2007h . This, along with 
the available surface pressure and feedbacks from HD 200775, sug- 
gests the six molecular cores resolved by our observations have the 
potential to collapse to form stars. 

4.2 Global motions 

Comparison of spectra of subregion / (Figure [4] (b)) with subre- 
gions // and /// (Figure [4] ( c) and (d)) evidently shows variation 
of line profile. Obviously, the line centers for the western subre- 
gion / are slightly blue-shifted with respect to the system velocity, 
while red-shifted for the eastern subregions // and ///. This velocity 
variation is also detected with 13 CO 7 = 2-1 line emission. The 
velocity centroid shows a tendency of increasing with ascending 
right ascension (J2000) for the six cores (see Tablejl). In Figure[8] 
the distribution of velocity centroid for the line of 13 CO 7 = 2 - 1 is 
presented and systematic velocity variation is confirmed. The value 
of the velocity gradient is estimated t o be about 0.94 km s~ ' pc~', 
which is consistent with the result in iGoodman et all jl993l) . This 
could be due to the strong stellar wind from HD 200775 which is 
the brightest source in this region. On the other hand, this velocity 
gradient may suggest a scenario of cloud rotation with a NW-SE 
direction, while the southeast potion moving away from us and the 
northwest portion to us. 

Feedbacks from stars with masses > 8 M e would profoundly 
affect the conditions of the natal clouds. The strong radiation pres- 
sure and stellar winds could destroy the clouds to constraint further 
star forming activities, otherwise they may help dense cores form in 
surroundings which will collapse to form new generation stars (see 
Sect l4.3t . In this work, HD 200775, a Herbig Be star with mass of 
~ 10M o , would generate violent winds blowing the whole molecu- 
lar cloud to globally expand. The powerful stellar winds or outflow 
from this young star blew away ambient material and excavate a 
cavity at the west where the density was relatively low. The dis- 
persed material accumulated around the cavity to form filamentary 
structures where core 1 and core 2 have formed. For the portion to 
the east of HD 200775, large amount gas and high density prevent 
it being blown about. However, the winds still possibly make the 
whole molecular to be globally expanding. 



Another observational evidence supporting a scenario of 
global expansion comes from the line profiles of 12 CO 7 = 3-2. 
As shown in Figure []](£>), red-skewed profiles dominates in spec- 
tra with two peaks. As mentioned in Sect l3.2l this is due to self- 
absorption of the molecular cloud with motions. Extensive red- 
skewed spectra can be observed in regions with expansion or in- 
clined outflows. In our observations, the distribution of the red- 
skewed lines is more consistent with the situation of expansion. 
Given the fact that an outflow always has lobes constrained within 
certain directions, red-dominated profiles originated from an in- 
clined outflow would be located in a region with relatively small 
opening angle and better collimation. In a scenario of an expand- 
ing molecular cloud with a hot bright star inside, the blue-shifted 
emission originate from the hemisphere closer to the observer while 
the red-shifted emission from the hemisphere farther away; for the 
excitation of the molecules is the highest near the hot star, the ob- 
server is always looking at the cooler side of the blue hemisphere 
and the hotter side of the red hemisphere; thus, the red emission 
should always be as strong as, or stronger than, the blue emis- 
sion (an analogue to the scenario of collapsing core proposed by 
IZhouetal.lfl993l) ). 

However, not all lines of 12 CO 7 = 3-2 emission resemble 
the red-skewed profiles. This ostensibly conflict with the expansion 
scenario. Blue-skewed spectra, rather than red-skewed ones, dom- 
inates at the northwest (see Figure [T]and|4j. This could be due to 
the infall motions in core 1 as aforementioned. The 12 CO 7 = 3-2 
spectra detected at the southwest show mono-peak with relatively 
large wing emissions (see Figure Q](fo)). This indicates that the gas 
column density is pretty low. Thus the line of l2 CO 7 = 3-2 here 
is optically thin while compared with the one in other places in 
this region. Another confusion comes from the complicated spec- 
tra in subregion /// as shown in Figure Q] where blue- and red- 
skewed line profiles are alternatively detected. Other mechanisms 
are necessary to destroy the homogeneity of 12 CO 7 = 3-2 line 
spectra. The blue-skewed profiles observed in the first column from 
right in Figure [4] (d) could be due to the backflow of stellar wind. 
To the southeast of the cloud, there are two 2MASS YSO candi- 
dates, 2MASS 15 (2MASS J2 1025484+68062 10) and 2MASS 16 
(2MASS J21025943+6806322, associated with IRAS 5). As shown 
in Figure|2]and|4](flJ, the cloud is as if compressed by the feedback 
from these infrared sources. With particularly large infrared exces- 
sive emission, 2MASS 15 could be in an early evolutionary stage 
with possible outflow activities which may contribute to the occur- 
rence of the blue-skewed profiles in the left three columns in Figure 
®(d). 

4.3 Triggered star formation 

Shown in Figure |9]is a close up view of the central part of the re- 
gion. Gray scaled background with an inverted color map presents 
4.5 ^m emission which mainly originate from H 2 (v = 0-0, 
s(9, 10, 11)) and CO(v = 1 - 0) in shocked regions or out flows 
dSmith et al.ll2006l ; lDavis et alj2007l ; lcvganowski et"ai]|2008l) . Ev- 
idently shocked structures are resolved at north and south of HD 
200775. It is believed that supersonic stellar winds and/or outflow 
originated from HD 200775 have been serving as the engine of 
these shocks. Relatively weak extended 4.5 /jm emission were de- 
tected in the east reaches. Nearby IRAS 4 (i.e., IRAS 21015+6757), 
an elongated structure still can be identified. This would be a trail 
feature of earlier shocks. For this region, we recommend a scenario 
of intermittent shocks generated by stellar winds from HD 200775. 
Overlaid contours on 4.5 jim emission in Figure [9] present 
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Figure 8. Distribution of line centroid of 13 CO 7 = 2-1. The overlayed contours are as the ones in Fig|2] 
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Figure 9. Close up view of the central part of the L 1 1 174 region. The gray 
scaled background with a inverted color map presents the 4.5 fim emission 
observed with PACS on Spitzer. The contours shows the velocity-integrated 
intensity of 13 CO 7 = 2-1. The crosses, filled triangles, and asterisks 
indicate the molecular cores, IRAS sources, and 2MASS YSO candidates 
respectively. The bright Herbig Be star HD 200775 is marked with a red star. 
The dotted ellipse encircles the 9 YSO candidates nearest to HD 200775. 
These potential YSOs are aligned in NW-SE direction. 



velocity-integrated intensity map of 13 CO 7 = 2-1 which roughly 
represents density distribution of gas. One can see large intensity 
gradients in sides toward HD 200775 for core 3, 4, and 5 which 
are located in sites closely surrounding HD 200775. These steep 
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Figure 10. 2MASS index as a function of the angle distance to HD 200775 
for the central 10 YSO candidates. 



increases of intensities from HD 200775 to the near cores could 
be due to the strong winds from the bright Herbig Be star. Fierce 
winds have been enhancing gas density around HD 200775 to help 
form cores which would collapse to form new generation stars. If 
this is the case, one would expect to see a sharp velocity gradient 
around the shell excavated by HD 200775. As one can see in Figure 
[8] the velocity gradient increases some around HD 200775, but not 
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sharply. This could be due to the poor spatial resolution of our ob- 
servation. The beam of 80" would smooth the velocity difference 
in 0.2 pc. 

Ten 2MASS YSO candidates (i.e., 2MASS 1-9 and 11) are 
marked with blue asterisks in Figure [9] They are labeled with num- 
bers with increasing distances from HD 200775. Note that the near- 
est nine ones are aligned in NW-SE direction consistent with the 
velocity gradient shown in Figure [8] These features are suggestive 
of a scenario of triggered star formation in this region. 

We define 2MASS index for a 2MASS source as, 



[J-H] 



l.8[H-Ks]- 0.1035 



(9) 



where [7 - H] and [H - Ks] are 2MASS color indices. 2MASS 
indices for all 17 YSO candidates were calculated and tabled in 
column 9 of Table [4] With an a = 1, a 2MASS point source 
will locate on the right dashed line in Figure [5] which is the 
dividi ng line between reddening normal field stars and T Tauri 
stars jBessell & Bred Il988l : iMever. Calvet. & Hillenbrand! 1 19971 ; 
iLi & Smith! 120051) . A younger YSO would possess a smaller 
2MASS index. 

We dotted the 10 2MASS YSO candidates nearest to HD 
200775(i.e., 2MASS 2-11, with exclusion of 2MASS 1 which is 
the 2MASS association of HD 200775) on a plot of 2MASS in- 
dex vs. angle distance to HD 200775 (see Figure |9J. A tendency 
of decreasing of a with increasing angle distance is presented. To 
quantitatively evaluate the existence of such inverse correlation, we 
performed least square linear fitting. Presented as a dashed line is 
the fitted result with all ten points. This is representative a linear 
equation of a = (-0.070 ±0.028)d+( 1.068 ±0.089), where d is the 
angle distance to HD 200775 from a 2MASS YSO candidate. This 
indicates that the more aged YSOs, with larger 2MASS index a, 
resides closer to HD 200775 and younger ones farther away. How- 
ever, a modest correlation coefficient of R 2 = 0.44 indicate that this 
tendency is not highly reliable. Shown as a solid line in Figure [9] 
another fitting without considering 2MASS 9 and 10, which pos- 
sess extremely small 2MASS indices, results in a linear equation of 
a = (-0.023 + 0.007)rf + (0.995 ± 0.021) and a larger correlation 
coefficient of R 2 = 0.62. Intriguingly, the tendency of descend- 
ing a (i.e., decreasing age) with increasing angle distance is better 
confirmed. This feature is reminiscent of sequential star formation 
which is always i nterpreted as evidence of t riggered star formation 
in the literature dElmegreen & Lad"alll977l ; IPreibisch & Zinneckerl 
120071) . 



5 SUMMARY 

We carried out mapping observations of LI 174 in lines of 13 CO 
7 = 2 - 1 and 12 CO 7 = 3 - 2 at 220.399 GHz and 345.796 GHz. 
Based on observed and archival data , the molecular conditions and 
star forming activities in LI 174 were extensively investigated. 

Six molecular cores are resolved by the observations in 13 CO 
7 = 2-1. Local thermodynamic equilibrium based calculations 
resulted in core masses of 5 to 31 M Q with an median value of 15 
M . Relatively large line width in all six cores indicates turbulent 
motions throughout the LI 174 region. This also leads to large virial 
masses ranging from 171.2 to 463 M , a factor of magnitude higher 
than the LTE masses. However, these cores are still possible to col- 
lapse to form new stars given the available external pressure and 
the affects from the bright Herbig Be star HD 200775. 

Conspicuously asymmetric line profiles of 13 CO 7 = 2-1 



spectra are observed in a large area in LI 174, which indicates 
bulk motions in this star forming region. Purely blue-skewed line 
features detected in core 1 suggest the existence of infall mo- 
tions. Prominent wing emissions imply potential outflows in LI 174 
which are not resolved because of the low resolution of our obser- 
vations. 

Extensive red-skewed line features are detected to the east of 
HD 200775. This, together with the velocity gradient along NW-SE 
direction, indicates global expansion in LI 174. Formed in the natal 
cloud, HD 200775 generated strong winds that blew away ambient 
material and excavated a cavity at the west where the density was 
relatively low. For the portion in the east, large amount of gas and 
high density prevent it from dispersing. And the consequence of the 
interaction between the winds and the gas is the global expansion 
of the reconstructed cloud. 

Seventeen YSO candidates were identified according to 
2MASS colour indices. For the ten 2MASS YSO candidates near- 
est to HD 200775, a tendency of decreasing age with increasing 
distance to HD 200775 is detected. This, along with large 13 CO in- 
tensity gradient of cores near HD 200775 and the shock features 
observed at 4.5 yum, suggests a scenario of star formation triggered 
by the Herbig Be star HD 200775. 
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